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Abstract 

A simple model based on the concept of resonant tunnelling clusters 
of lattice defects is used to explain the low temperature anomalies of hep 
''He crystals (mass decoupling from a torsional oscillator, shear modu- 
lus anomaly, dissipation peaks, heat capacity peak). Mass decoupling is 
a result of an internal Josephson effect: mass supercurrent inside phase 
coherent tunnelling clusters. Quantitative results are in reasonable agree- 
ment with experiments. 
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1 Introduction 

Helium crystals are the most pronounced examples of the quantum crystals 
in which quantum tunnelling of particles results in many unusual phenomena. 
Delocalization of impurities, vacancies, dislocations, and surface defects leads 
to quantum diffusion, vacancy-induced mobility of impurities and ions, internal 
friction anomalies, and weakly decaying crystallization waves, respectively (see 
the review paper |T|). 

Motivated by theoretical predictions l2|(3','4' of the superfluidity of solid ^He, 
Kim and Chan (KC) performed experiments [5, 6^ similar to the Andronikashvili 
experiment [71 18] in which the decoupling of the superfluid fraction of liquid Hell 
from a torsional oscillator (TO) was discovered. The remarkable observation of 
KC was a similar decoupling of a part of solid helium from TO below 0.2 K 
which was interpreted as the superfluidity of a solid. 

However, unlike the superfluid transition, the onset of the mass decoupling 
is broad and is accompanied by a dissipation peak. Day and Beamish 9 mea- 
sured the elastic shear modulus, and observed a similar behavior (stiffening and 
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dissipation) in the same temperature range. Near the decoupling onset temper- 
ature, a broad heat capacity peak was observed [10], but no pressure- induced 
superflow through the sohd was found [TTl [T^] . The magnitudes of anomaUes 
depend strongly on the way the solid was prepared [T?. All these anomalies 
seem to be absent in perfect crystals. 

In this paper (see also earlier letter [2]), a simple model is proposed to 
explain the low-temperature anomalies of imperfect *He crystals. As in papers 
[TB] . the concept of tunnelling two-level systems (TLSs) in solids [T7] is 
used. In [TSl [H] TLSs were considered in highly disordered (glassy) samples. 
In this case, the parameters of TLSs are uniformly distributed according to 
the original tunnelling model |18| . However, in the most experiments, solid 
^He samples were grown by the blocked capillary technique, but they consisted 
of poly-crystals. TLSs in crystals are degenerate (or resonant) TLSs in which 
the bare (with no tunnelling) energy difference of two localized states is zero 
(see [T7] and below where a simple example is presented). The physical reason 
of the degeneracy is the crystal symmetry. Two localized states transform to 
each other under a crystal symmetry transformation. (Otherwise, there is no 
reason for energy difference to be small in crystals) . The degenerate TLS is only 
the simplest case. Tunnelling systems may generally consist of more than two 
localized configurations |17) . 

In our simple model, we suppose that the main contribution to all anomalies 
is introduced by degenerate TLSs of a certain (or crystallographically equiva- 
lent) structure. In this case, the characteristic parameters, including tunneling 
amplitudes, are the same for all TLSs. 

The key point is a peculiar quantum phenomenon of momentum deficit for 
TLSs in moving solids [HI [16]. In a solid with local velocity v, the momentum 
of a TLS can under certain conditions (see below) be equal to m*v, where m* is 
the effective mass which is different from the contribution m of the TLS to the 
total mass. Generally, m* depends on the frequency and amplitude of the local 
velocity oscillations, the difference m* — m being always negative. Therefore, 
the TLS has nonzero internal momentum (m* — m)v directed opposite to the 
velocity of the solid. In classical physics, a system of particles moving in a 
restricted spatial region, has zero internal momentum up to the frequencies on 
the order of the characteristic frequency of particle oscillations. For TLSs, this 
frequency is the tunneling frequency. We show below that under conditions of 
KG experiments, TLSs in solid helium have the frequency independent internal 
momentum down to the frequencies six orders of magnitude lower than the 
tunneling frequency. The internal momentum disappears at frequencies below 
the reciprocal transverse relaxation time (the TLS phase memory time). 

2 TLS in imperfect crystals 

Simple example of a resonant TLS is a four- vacancy cluster in a hep ^He crystal 
(see Fig.l). Vacancies are located in the apexes of a tetrahedron. Three of them 
are disposed in the symmetry plane, which is perpendicular to the c-axis of the 
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Figure 1: Hypothetic arrangement of vacancies in a tetrahedronal cluster with 
possible tunneling "through" the symmetry plane. 

crystal. The fourth vacancy occupies one of the two positions which transform 
to each other by the reflection in the plane. According to numerical calculations 
|19j . four- vacancy clusters in solid helium are either unbound or are bound too 
weakly at the temperatures of KC-experiments. However, highly metastable 
growth-introduced vacancy clusters can persist in crystals for a long time. The 
same can be true for clusters of a different metastable phase in hcp-crystals 
[20] . Large growth-introduced vacancy clusters are probably the faceted liquid 
bubbles studied experimentally in '^He crystals |21) . 

The Hamiltonian of a TLS in a stationary crystal is Hq = + Aai , where 
(7a, {a — 1, 2, 3), are the Pauli matrices, A is the tunneling amplitude, and 
{£, > 0) are the bare energies of two localized states. Here, the TLS is supposed 
to be nondegenerate, because we consider below a shear deformation that breaks 
the symmetry. 

3 TLS in rotating crystals 

The Hamiltonian of the TLS in a crystal moving with a local velocity v, accord- 
ing to Galilean transformation is Hi = iJo + pv, where p is the TLS momentum 
operator. We suppose that the TLS tunneling is accompanied by the displace- 
ment of a mass to by a vector a. The operator part of the TLS center of 
gravity coordinate can be written as r = — (T3a/2. The momentum operator is 
p = mr = {im/h)[Hor] — —{mA/h)aa2 [TS l fTB ] . 

The internal momentum in the TLS state ipi |1) -I- "02 |2) is proportional to 
((T2) oc sin((pi — ip2)- Here, |1) and |2) are the localized states, and ipi and (p2 
are the phases of t/'i and 1^2, respectively. Therefore, the internal momentum 
is caused by a kind of internal Josephson effect: the mass supercurrent inside 
phase coherent TLSs. 
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We suppose that the velocity is a result of an axisymmetric container rota- 
tion. Otherwise, additional terms should be added to the Hamiltonian to take 
the macroscopic displacement of the container walls into account (see [H], §11). 
Because the size of the TLS is assumed to be much smaller than the length scale 
of the rotating container, we can use the following expressions for the velocity 
and the TLS angular momentum: v = x R, and M = (R x p)j,, where is 
the angular velocity and R is the coordinate of the TLS (R ± il) with respect 
to the origin located at the rotation axis which is parallel to the z-axis. We 
have 

Hi=Ho + MQ, (1) 

where M = fj,a2, ^ = —{mA/h)aRcosd, 9 is the angle between a and v. 
After canonical transformation by the unitary operator 

C/ = (2e(£ + 0)-'/'(£ + e-»A'T2), (2) 

where e = (A^ +^^)^^^, we obtain the Hamiltonian in the more convenient form 

H = UHiU+ = -e(j3+ nVL(j2. (3) 

This Hamiltonian can be written in the form 

H = -haaa, (4) 

where ha is the "field" with components hi = 0, /12 = — M^^, and hs = e. 

The TLS density matrix w is generally determined by the real polarization 
vector Sa- 

= (l + s„<Ta)/2. (5) 

We have 

(ctq) = Tr (wo-q) = Sa- (6) 
From the equation for the density matrix 

w^{t/h)[w,H], (7) 

we obtain the dynamic equation for a free TLS (without dissipation): 

hsa = eapjhps^, (8) 

where Cap^y is the Levi-Civita tensor. 

The adiabatic theorem (see [23], chap II, §5c) takes place as a consequence of 
(8): along with the modulus s = \sa \ of polarization, the angle between the field 
ha and Sq is the integral of motion. The process is adiabatic if the time scale of 
the field variation is much longer than h/\ha\. The last condition is fulfilled with 
a significant safety margin: the characteristic frequency of the field variation in 
experiments mentioned above, is on the order of 1 kHz, but the magnitude of 
the field \ha\ is on the order of A ~ 0.1 K. 

Until the rotation (and deformation) is applied, the polarization is directed 
along the field, and the absolute value of the equilibrium polarization is s = sqi 
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where sq = tanhe/T. According to the adiabatic theorem, we have for the free 
TLS: 

Si=Q, s,^--ns, .3^_,s^.__^_-^., (9) 

1 /2 

where i? = (e^ + ^^fi^) is the field modulus. 

We suppose that equations describing relaxation processes in a stationary 
solid have the standard form: 

S2 . S3 -So 
S2 = , S3 = , (10) 

where ti and T2 are the longitudinal and transverse relaxation times, respec- 
tively. Additional terms to the time derivatives (10) taking the rotation of the 
solid into account, can be determined from (9). Finally, we obtain the following 
dynamic equations for the TLS: 

d f \ S2 d( L?VP_\ S3 - So 

Following the work of Burin et al. [23], we will assume that both reciprocal 
relaxation times depend linearly on temperature: 

rf'-XiT, T^^^X2T, (12) 

where xi ^-^d X2 are constants. 



4 Torsional oscillations 

Consider torsional oscillations with a small amplitude VL(t) oc exp(— iwt). From 
the first of the equations (11), we obtain the mean value of the TLS angular 
momentum (Af ) = /iS2 = /(w)f^, where 



11^ iujT2 - uj^rj A 



is the TLS rotational inertia. The real part of / determines the TLS 

mass deficit 5m{uj) = m — m* = —I'{lo)IE?. The crystal fraction /, decoupled 
from TO, is / = n5m/ p, where p is the crystal density, and n is the number of 
resonant TLSs per unit volume. We have 

^ , , f ma cos uj'^t^ , , A ,^ 

M-)=(^^— j ^;;^Atanh-. (14) 

The mass deficit of a free TLS (a;r2 ^ 1) is frequency independent and coincides 
with the result of [TBI. 

For uniform (w = 0) rotation Sm = 0, in accordance with |15l 116) . The 
intermediate "quasi-equilibrium" region which was considered in [151 IS], is 
absent in our simple model. 
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Korshunov [^S] calculated the TLS mass deficit for translational oscillations 
of a closed container filled with solid helium. His value of Sm is proportional 
to the extremely small factor (w/A)'^. In fact, the Korshunov result is an ex- 
planation of the absence of TO anomalies in the blocked annulus experiments 

The contribution (E) — (flM) of the TLS to the energy dissipation is deter- 
mined by the imaginary part /"(w) oil (see [Hj, §123): {E) ^ (a;/2)/"(a;) 
Under experimental conditions [27|, the dependence of the reciprocal quality 
factor of TO on temperature and frequency, is determined by the expression 

(15) 

The temperature dependence of (15) for two different frequencies cj are shown 
at the bottom of Fig. 2. Note that according to (15), these curves intersect at 
some point. 

The corresponding formula for the decoupled fraction /, according to (14), 

is 

^°^:;5^tanh-. (16) 

The temperature dependencies of (16) for the same two frequencies are shown 
at the top of Fig. 2. 



5 Nonlinear torsional oscillations 

Consider nonlinear torsional oscillations at low temperatures (free TLSs). The 
oscillations are described by the equation 

& + uj^a^ -{l/Io)M{t) (17) 

for the rotation angle a {a = fl). Here, Iq and ujq are the rotational inertia and 
the TO fundamental frequency in the absence of TLSs, respectively. M{t) is 
the TLS angular momentum. According to (1) and (9) we have 

i?A wcos^ 6 

" (1 + cos2 e{v/vc)^y/^ ' ^ ' 

where Vc — h/(ma) is a characteristic (critical) velocity. 

In the absence of TLSs, the oscillations are harmonic: a{t) ~ ao cos ip, where 
ao is the amplitude and ip = luqI. The time dependence of M is determined by 
(18) with V = vosinip, vq = — agWo-R- Being the periodic function of ip, M{t) 
can be represented by its Fourier series. It is known |28| that a small correction 
to the fundamental frequency is determined by the resonant part (proportional 
to cos ip) of the right hand side of (jl7p . Resonant contribution to the momentum 
M,:es{ip) is proportional to simp 

1 f" 

Mi.cs(<p) = sin(^- / M{ip')smp'dip' . (19) 
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Figure 2: Temperature dependencies of decoupled fraction (|16p and dissipa- 
tion (|15p (A — 51 mK, X2 = 5.9kHz/K) for two frequencies fitting experimental 
data [H]. 
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According to (18) and (19), the TLS contribution to the fundamental frequency 
and, therefore, to the value of the decoupled fraction f{vo) is determined by the 
expression 

/IT ■ 2 

(l + cos^^(t;o/Uc) sm v') ' 

We consider a narrow annulus where R » const for all TLSs. For poly-crystals, 
the expression (20) should be averaged over all directions of a at a given v. We 
have 

f{v) = const cos^ edcos0 f ''f ^f'^ . , (21) 

J -I J-T, (1 + {v/VcY cos^ ^sm^ <^)V^ 

Here and below we omit the index in the velocity amplitude vq. Finally, we 
obtain 

/(O) ~ TT io io (1 + {v/v.?x^ sin^ ^)i/2 • 
This dependence is plotted in Fig. 3. According to (14), the mass deficit 5m for 
a single TLS is on the order of A/v'^. Assuming A ~ 0.1 K and Vc ~ 10 /im/s, 
we obtain 6m ~ 10~^^ g. The experimental value of the decoupled fraction / = 
n5m/p ^ 10~^ corresponds to the quite low TLS concentration n ^ lO^'cm^'^. 
The characteristic size L of the defect clusters can be estimated from the value 
of the critical velocity. Assuming that ma = pL^, we have L ~ 10 nm. 



6 The shear modulus anomaly 

Let C = dux/dy be a shear strain (u^ is the displacement along the a;-axis). 

According to (3) in the absence of rotation the derivative of the TLS Hamiltonian 
with respect to the strain is dH/dQ = —azde/dC,. The TLS contribution to the 
stress a is determined by the mean value of this derivative 

a = n< dH/dC >= -nssde/dC,. (23) 

We consider small oscillations of the strain C,{t) oc exp(— iwt). The oscillat- 
ing part of the TLS polarization Sg and that of the instantaneous equilibrium 
polarization s'q satisfy, according to (10), the equation iujs'^ = (sg — s'q)/ti, 
where 

s^ = C^tanh|; = 1 . (24) 

Here, {de/dQo is the derivative de/dC, at C = 0. We have 



i C (de\ 

1 - iojTi Tcosh^(A/T) V^C/o 
The oscillating part of the stress is 

-' = -4cfS') tanh^+fg^) s'A. (26) 
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Figure 4: Experimentally observed (30j temperature dependence of shear mod- 
ulus at high frequency and a fit for A = 0.08 K. 
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Figure 5: Temperature dependence of heat capacity ((^ (A — 0.11 K) and 
experimental data [lOj . 

The TLS contribution to the complex shear modulus G{uj) defined by the ex- 
pression a' — G{uj)(, is 



G(w) 



2^ 



, A 1 1 fde\ , 

tanh 5 — y (27) 

T l~iLjnTcosh\A/T) \dCj ol 



The second derivative d^e/dC,"^ should be negative. Otherwise, the real part of 
G, contrary to experiments [9] , would be negative. 

The dissipation of the shear oscillations is determined by the imaginary 
part of G. Its temperature dependence is characterized by the presence of a 
characteristic peak, in accordance with experiments ^ . 

It was noted [9] that G' and / change very similarly with temperature. This 
is consistent with our results: due to (16) and (27), in high frequency region 
ujT ^ 1, we have G oc / oc tanh(A/T). This dependence for G is plotted in 
Fig.4. 

7 Thermodynamics of resonant TLSs 

The contribution of resonant TLSs to the free energy of the unit volume of a 
crystal is 

A' 



F = -nT log ( 2 cosh — ) . (28) 
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The heat capacity is determined by the well-known formula: 



( 



A 



) 



2 



C = -T 



QT2 



n 



T 



cosh ^ — 
T 



(29) 



This dependence is plotted in Fig. 5. 

The contribution of TLSs to pressure is 




(30) 



The P{T) dependence at a fixed volume was measured by Grigoriev at al. [3T| . 
The glassy contribution (a T^) was found. The term (30) was absent. The 
possible reason is that the derivative dA/dp is anomalously small. Indeed, 
this derivative determines in our model the pressure dependence of the onset 
temperature of the mass decoupling. KG [29] observed no apparent change of 
the onset temperature with pressure. 

8 Summary 

Low-temperature anomalies of hep ^He crystals are caused by growth-introduced 
resonant nanoscale clusters of lattice defects. The quantum phase coherence and 
tunnelling of the clusters are crucial. The important problem is to identify the 
structure of the clusters. The most probable candidates seem to be the faceted 
liquid bubbles. 
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